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DESIGN PROCEDURES 
FOR 

FLUTTER-FREE SURFACE PANELS 


Robert M. Laurenson and J. I. McPherson 
McDonnell Douglas Astronautics Company - East 
St. Louis, Missouri 

SUMMARY 

An approach for the design of lightweight external surface panel con- 
figurations to preclude panel flutter has been developed. Design procedures 
were developed for flat orthotropic panels under the interacting influence 
of parameters such as support flexibility, inplane loads, pressure differen- 
tial, and flow angularity. The relationships required to define these 
design procedures were based on theoretical panel flutter analyses. Where 
possible, the design procedures ivere verified through comparison with avail- 
able experimental panel flutter data. 




1.0 INTRODUCTION 


Future high velocity cruise aircraft and reusable entry vehicles must 
cope with environments which place severe thermal and structural demands on 
the external surfaces. Prevention of panel flutter in these external sur- 
faces is one of many problems which designers must face. Traditionally, 
external panels have been designed by strength, buckling, and fatigue con- 
siderations and then checked for panel flutter performance. Such a procedure 
often leads to "flutter fixes" and additional weight for panel flutter pre- 
vention. 

To permit consideration of panel flutter early in the design process, 
design procedures have been developed in a format useful to designers without 
prior panel flutter experience. The major emphasis is placed on development 
of flutter design procedures and their application to metallic panels with 
orthotropic stiffness characteristics and a variety of edge support conditions 
from nearly free to fully clamped. The panels may also be subjected to 
various loading and flow conditions. Additionally, results are also presented 
for surface panels attached to primary structure by an elastic coupling medium. 
For either panel configuration it should be emphasized that although these 
design procedures are useful in preliminary design, detailed flutter analyses 
and/or tests may be necessary to verify the final design of critical panel 
configurations. 

A wealth of theoretical panel flutter investigations and wind tunnel 
data exists in the literature. The theoretical approaches were employed to 
establish the trends and relationships required to define these panel flutter 
design procedures. Experimental flutter data were used to verify the de- 
sign criteria which are presented in the form of flutter-free panel design 
boundaries. The flat, orthotropic panel configuration to which this approach 
has application is illustrated in figure 1. The influence of such interacting 
parameters as flow angularity, support flexibility, inplane loading, ortho- 
tropic panel stiffness characteristics, damping, and static pressure differen- 
tial on the flutter of external surface panels is included in these procedures. 
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FIGURE 1 SURFACE PANEL CONFIGURATION 

Interest in the establishment of panel flutter design procedures has 
been high for a number of years. An initial attempt to define such procedures 
is given in reference 1-. In that document, a panel flutter design boundary 
is defined which envelopes the available experimental data from both unstiff- 
ened and corrugation stiffened panels. 

Modifications and refinements to this design approach have continued 
(refs. 2 and 3). In 1964, and later in a revised form in 1972, the National 
Aeronautics and Space Admir.istration issued a formal design monograph (ref. 4) 
to be used as a guide in the formulation of design requirements and specifi- 
cations dealing with panel flutter. Reference 4 does not provide detailed 
procedures to be followed for the design of flutter free surface panels but 
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rather the philosophy to be followed during the design process. Reference 5, 
published in 1968, gives simplified criteria in graphical form for many of the 
parameters important for panel flutter design; however, effects of panel 
orthotropy, edge support flexibility, and damping are not included. The panel 
flutter design procedures presented in this report are extensions of the 
previously mentioned design approaches. Additional parameters such as panel 
orthotropy, edge support flexibility, flow angularity, and damping, as under- 
stood in the current state of the art, are included in a graphical format 
which should be useful to designers without prior panel flutter experience. 
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2.0 SYMBOLS 


Definitions of the major symbols required for the application of the panel 
flutter design procedures are given below. Symbols not having general usage 
throughout this report are defined as they are introduced. 


a, b 
D 

^f 

^ 2 * ^12 
E 


c 

f(M) 

FP 


Panel length in x and y directions, respectively; m 
Isotropic panel stiffness parameter; N-m 
Damping factor - equation (41) 

Orthotropic panel stiffness parameters; N-m 

p 

Young's modulus of panel material; N/m 
Panel edge rotational restraint coefficient 
Mach number correction factor - figure 7 
Flutter Parameter 


g 

9t 

G 

GP 


h 


h 


eq 




* 


K 


If * 


M 


Structural damping coefficient 
Total damping coefficient 

2 

Shear modulus of panel material; N/m 
Geometry Parameter 
Isotropic panel thickness; m 

Equivalent isotropic panel thickness - equation (39); m 

Nondimensional leading and trailing edge support stiffness 
parameter 

Nondimensional stiffness parameter IT. ^ modified for presence of 
flow angularity - equation (29) 

Nondimensional stream edge support stiffness parameter 

Nondimensional stiffness parameter 1c modified for presence of 
flow angularity - equation (28) 

Mach number 

Critical buckling load; N/m 
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Inplane load in x direction; N/m 

Nondimensional pressure differential parameter - equation (38) 
Pqp Buckling load ratio - equation (12) 

q Free stream dynamic pressure; Pa 

Q Pressure differential correction factor - figure 6 

r 

S|^y Leading and trailing edge support parameter - equation (22) 

* 

Support parameter modified for presence of flow angularity 
- equation (30) 

t T ime 

a Equation (4) 

Y Panel mass density per area 

Ag Differential damping coefficient - equation (42) 

Ap Static pressure differential; Pa 

X Flow parameter - equation (6) 

V Poisson's ratio of panel material 


The International System of Units 


port. The appropriate quantities are 

Quantity 

Unit 

length 

meter 

force 

newtpn 

pressure 

pascal 

Abbreviations for 

the following 

of units in this report 

• 

Prefix 

Multiplying 

Factor 

mill i 

10"^ 

kilo 

10^ 

mega 

10® 

giga 

10^ 


(SI) has been used throughout this re- 
ined as follows: 

SI Symbol 
m 
N 

Pa 

ixes have been employed for multiples 

Abbreviation 

m 
k 
M 
G 
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3.0 FLUTTER-FREE PANEL DESIGN BOUNDARIES 


To allow designers without prior panel flutter experience to consider 
panel flutter early in the design process, flutter-free panel design boundaries 
have been developed in terms of a nondimensional panel geometry parameter GP 
and a nondimensional flutter parameter FP. Effects of interacting parameters 
such as panel size, panel stiffness characteristics and panel support condi- 
tions are included in the geometry parameter. The flutter parameter contains 
the effects of free-stream dynamic pressure and Mach number. 


The initial definition of the geometry and flutter parameters and their 
relationship was based on a two mode solution for flutter of simply supported 
orthotropic panels. The following assumed panel deflection satisfies simply 
supported boundary conditions 

w(x,y,t) = W^.| (t)Sin Sin -^ + Wgi (t)Sin Sin ^ (1 ) 


The basic partial differential equation for flutter of an orthotropic 
panel (ref. 6) is of the form 


D — + 

w 


2D 







( 2 ) 


where the aerodynamic loading has been assumed to be given by two-dimensional 
static aerodynamics. Employing the assumed panel deflection of equation (1), 
a Galerkin solution to the above yields the following set of ordinary 
differential equations: 


(?/ [ “i * (l/ " “2 (l) " ] “n - 3 ^ “21 * ^ "n 


= 0 


(3) 


. 16q 

^21 ^ 3fl^ ”11 


w„ + Y W21 = 0 
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Assuming a solution to equation (3) of the form 


and solving the corresponding characteristic equation yields the following 
two mode flutter solution 


, - 3tt^ 

^ ■ TT 


(t) 


(5) 


In equation (5), the flow parameter x is expressed as 

3 


X = 


- 2qar 


d{7(m) 


( 6 ) 


The geometry and flutter parameters are arbitrarily defined as 



Combining equations (7) and (8) with equation (5) yields the following 
relationship between GP and FP for simply supported panels 


= 0-0365 

5+2 (GP)^ 


(7) 

( 8 ) 


(9) 


A second approximate flutter solution, referred to as the preflutter 
solution, has been presented in several documents such as reference 7. This 
approximate solution is very accurate for simply supported panels for values 
of GP > .1 and results in the following relationship between GP and FP 


FP 


0.0593 

(5 + GP^) V 4 + 2 GP^ 


( 10 ) 
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Results from the two mode solution and the preflutter solution are 
presented in figure 2 where GP is shown as a function of FP. Both solutions 
have the same trends; however, the two mode solution lacks sufficient accuracy 
for design purposes and the preflutter solution is not applicable to panels 
with edge supports which approach a fully clamped condition. 

Since the approximate solutions were insufficient for design purposes, 
reliance was placed on more detailed flutter analyses of panels with boundary 
condition combinations which varied from nearly free to fully clamped to 
obtain the final flutter- free panel design boundaries shown in figure 3. These 
analyses included both closed form solutions and well converged modal solutions 
so that accurate results were used to develop the boundaries. 

As the effects of different parameters such as inplane loading, flow 
angularity, and support flexibility were considered it became necessary to re- 
define the geometry and flutter parameters to maintain the relationship between 
these two parameters given by the flutter- free panel design boundaries shown 
in figure 3. The required expressions for GP and FP were determined by combined 
analysis and curve fitting of experimental data. This approach was followed 
to duplicate the trends obtained in detailed theoretical and experimental 
panel flutter investigations reported in the literature. 

As an example, the two mode flutter solution was extended to include 
the influence of inplane loads. This analysis resulted in a definition for 
GP to be used when inplane loads are present. To be specific, it can be 
shown that the two mode solution, including the influence of an inplane load 
in the direction of flow, N.| , takes the form 

( 11 ) 
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GEOMETRY PARAMETER - GP 



FLUTTER PARAMETER - FP 


FIGURE 2 FLUTTER-FREE PANEL DESIGN BOUNDARIES FROM APPROXIMATE SOLUTIONS 





GEOMETRY PARAMETER 



In equation (11) is the ratio of the inplane load N-j to the buckling load 
Thus, 


^CR “ 


( 12 ) 


Comparing equations (5) and (11) reveals that the second term of 
equation (11) indicates a new definition for GP of the form 





(13) 


As will be discussed in Section 4.3, equation (13) is in fact the form for GP 
that was selected to account for the presence of inplane loads. For a panel 
with inplane loads, modifications to the FP relationship of equation (8) were 
then made by correlating flutter predictions with available experimental data. 


For ease of computer application and to assist the designer in making 
accurate determinations of FP once GP has been calculated, analytical expres- 
sions of the design boundaries were determined by curve fitting techniques. 
These expressions take the form 


FP 


5+2 (GP)^ + 0.18 (GP)^ 


(14) 


where is the panel edge rotational restraint coefficient defined as 

E = 0.016 
c 


(15) 


for a panel with zero edge slope and 


= 0.027 (16) 

for a panel with zero edge moment. Thus for a given value of GP, the corres- 
ponding magnitude of FP may be obtained from equation (14) with the appropriate 
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The relationship of equation (14) is valid over the range of GP’s be- 
tween slightly more than zero (0.1) and 5.0. For values of GP lower than 0.1, 
the coefficient is defined as 

« 0.0157 (17) 

for a panel with zero edge slope and 

» 0.0292 (18) 

for a panel with zero edge moment. 
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4,0 DESIGN BOUNDARY APPLICATION 


The steps required in the application of the flutter-free panel design 
boundary during preliminary design of lightweight external surface panels are 
detailed in this section. Parameters which are included in this design 
approach as they influence panel flutter are: 

0 Panel aspect ratio 
0 Panel orthotropic properties 
0 Support flexibility 
0 Inplane loads 
0 Static pressure differential 
0 Flow angularity 
0 Structural damping 

Verification of these design procedures is addressed in Section 6.0. 


In the application of the flutter-free panel design boundary, the basic 
approach involves the determination of the magnitude of GP for the panel 
configuration of interest. The design boundary (fig. 3) is then entered 
with this value and the corresponding magnitude of FP obtained. With this 
value of the Flutter Parameter, the flutter critical flight condition may 
be defined. The forms taken by the parameters FP and GP depend on the 
particular situation under consideration. Definition of these parameters 
and their application in conjunction with the design boundary is detailed 
in the following sections. 


The steps required for application of the panel flutter design procedure 
during the design of lightweight surface panels are summarized in figure 4. 

The basic geometry data required are the panel length a, width b, and ortho- 
tropic stiffness properties Dp D 2 , and D^ 2 * ^ number of techniques for the 

determination of these panel stiffness parameters are presented in Appendix A. 
For an isotropic panel these three panel stiffness properties are equal and 
are expressed as 

h3r 

D = (19) 


12(l-v^) 
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FIGURE 4 PANEL DESIGN PROCEDURE 








When the influence of deflect! onal support flexibility is to be in- 
cluded in the design the nondimensional stream edge stiffness parameter 
and/or leading and trailing edge parameter Kj^j are required. These quantities 
are defined for a number of support conditions in Appendix B. 

Several options available for use with the flutter-free panel design 
boundary are indicated in figure 4. This is to aid the 'designer in select- 
ing the correct procedures to be followed when applying these panel flutter 
design criteria. For each option, reference is made in the figure to the 
report section detailing the application of the design boundary to a partic- 
ular panel configuration. A number of examples detailing the use of these 
panel design procedures are discussed in Section 5.0. 

The procedures to be followed in defining the magnitude of the geometry 
parameter for the panel configuration of interest are presented in the follow- 
ing sections. With this quantity, the flutter parameter is obtained from the 
flutter-free panel design boundary. This value of FP is then used to obtain 
the critical dynamic pressure parameter q/f(M), where q is the free stream 
dynamic pressure and f(M) is the Mach number correction factor. 

Modifications to the dynamic pressure parameter for a static pressure 
difference across the panel are covered in Section 4.4. As discussed in 
Section 4.5, the quantity q/f(M) may require modification to account for the 
presence of significant structural damping in the system. With this modified 
dynamic pressure parameter, the flutter critical dynamic pressure is obtained 
with the procedures detailed in Section 4.6. 

4.1 Basic Panel - For an orthotropic panel with no inplane load and 
zero flow angularity, the geometry and flutter parameters are defined as 


and 



C20) 


19 


FP = S 


f(M) 

*-T "7^3 
q a 


C21) 


In equation (20), is the nondimensional stream edge deflectional support 
stiffness and in equation (21) is the leading and trailing edge support 
parameter expressed as 


S 


LT 



( 


■2K 


LT 


(a/b)^ ( 0 ^ 2 /^^) 


) 


( 22 ) 


where Kj^j is the nondimensional stiffness parameter associated with the 
leading and trailing edge support flexibility. Additionally, the quantity 
C required in equation (20) is defined by the following expression 


V 02 


In this application the concept of support flexibility refers to a panel 
edge which is not completely restrained against deflections normal to the panel. 
Thus, in this case the panel boundary condition is less rigid than either a 
simply supported or clamped edge. Appendix B contains a number of procedures 
which may be employed to determine the magnitudes of the nondimensional edge 
support parameters and For panel configurations with unequal flexible 

supports on opposite edges an average value of the nondimensional edge support 
parameters 1<2 and l<yy should be used except for values of l<yy less than about 
five. At values of l<yy less than five a flutter analysis should be conducted 
since theory (ref. 8) indicates large reductions in flutter margins can occur 

at low values of l<yy for unequal supports at the panel leading and trailing 
edges . 

The theoretical trends presented in reference 9, of flutter critical 
flow conditions versus leading and trailing edge support stiffness, were 
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employed in obtaining the equation (22) definition of The modifying 

term present in equation (20), which accounts for the stream edge 

support stiffness, was obtained by matching the experimental trends presented 
in reference 10. 

As the stream edge deflectional support stiffness IL approaches 

p ^ 

infinity the quantity C /K^ approaches zero and equation (20) becomes 


Likewise, for a approaching infinity, the magnitude of approaches 
one and equation (21) may be written as 


FP 


f(M) 
q a^ 


(25) 


With the magnitude of GP defined with either equation (20) or (24), the 
corresponding value for FP is obtained from the design boundaries in figure 3* 
The dynamic pressure parameter q/f(M) may be determined as appropriate, from 
equation (21) or (25). Potential modifications to q/f(M) for static pressure 
differential are presented in Section 4.4 and for structural damping are 
covered in Section 4.5. The evaluation of f(M) allowing definition of the 
flutter critical dynamic pressure is discussed in Section 4.6. 


4.2 F low Angularity - For flow at an angle, and with edge flexibility 
included, the geometry and flutter parameters are defined as 



D^2/Di 


1 + C^/K 


2 * 
S 


Cos^ A 




Sin^ A 


(26) 


*Choice of the design boundary is dependent on the degree of rotational re- 
straint along the panel edges. Most panels are probably adequately repre- 
sented by the zero edge slope boundary; however, a more conservative design 
results from the zero moment boundary. 
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and 


FP 


f(M) 


qa 


2 

Cos A 



(27) 


In equations (26) and (27), A is the magnitude of flow angularity shown in 
figure 1 and and are edge support flexibility parameters modified 
to account for the flow angularity. 


The modified nondimensional edge support stiffness parameters including 
the influence of flow angularity are given as 


K$ = 


1 


Cos^ A Sin^A 


(28) 


and 


_ * 


Ks " Klt 


Sin^ A Cos^ A 


(29) 


+ K 


LT 


The modified leading and trailing edge support parameter is of the form 

/ 71 \V2 / \ 

V^\t''\t/ \(a/b)^ (D^2/°1^/ 




+ e 


+ e 


-2^'^LT (30) 

Sin 2A 


The expressions for GP and FP given by equations (26) and (27) were 
chosen to assure correct results for a flow angle of 90 degrees. That is, 
for a 90 degree flow angle, the a and b panel dimensions and the panel stiff- 
ness properties are properly interchanged. The combination of flow angular- 
ity and panel edge support flexibility has been expressed in terms of effective 
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edge support stiffness parallel and perpendicular to the flow. Thus, the 
modified edge support stiffness parameters of equations (28) and (29) are 
based on the assumption that the leading and trailing edge springs act in 
parallel with the stream edge springs. The last term of equation (30) is 
included to match the experimental trends for panels with flow angularity 
such as presented in reference 11. 

For panel configurations which include both edge support flexibility 
and flow angularity, the application of equations (28) through (30) in 
combination with the relations presented for FP and GP can result in inaccu- 
rate flutter predictions. This occurs when the second term in the expression 
for [eq. (22)] is not small when compared to the magnitude of the first 
term. To overcome this difficulty, an equivalent nondimensional leading and 
trailing edge support stiffness was defined as 




) 


+ 36 


-6 


(31) 


where is obtained from equation (22) for the particular value of 
The magnitudes of GP and FP are then obtained from equations (26) and (27) 
where and are based on the quantities and 


For configurations where the edge support stiffness parameters and 
approach infinity (total edge deflection restraint), equations (26) and 
(27) become 


and 


= tl D 


12 


Cos A + 


(.y s,. 


FP = 


f(M) 
q a^ 


Cos A + 




(32) 


(33) 
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It should be noted that the definitions for the geometry and flutter 
parameters and the edge support stiffness parameters are based on data at flow 
angles of zero and in the range from about 15 to 90 degrees. Analysis (ref. 6) 
indicates that slight variations in flow angle in the region of zero to 15 
degrees can have large effects on flutter conditions. Due to lack of experi- 
mental flutter data for flow angles of zero to 15 degrees some uncertainty 
exists about the validity of the current design approach in this region. When 
effects of flow angularity are to be included a conservative design can be 
obtained if the panel is considered for both zero and 90 degree flow and the 
design based on the worst case. 

With the magnitude of FP obtained for the appropriate deflection edge 
support condition, the dynamic pressure parameter q/f(M) can be evaluated. 

The effects of static pressure differential, Section 4.4; structural damping. 
Section 4.5; and Mach number, Section 4.6, should be examined and included in 
the determination of the flutter critical dynamic pressure when significant. 


4.3 Inplane Loads - The flutter-free panel design boundary has appli- 
cation to panels subjected to inplane loading if the flow is parallel to the 
X axis. In addition, the panel stream edges may have either total deflec- 
tional restraint or some edge support flexibility as defined by The 
leading and trailing edges are assumed to have complete deflectional restraint. 
Figure 5 illustrates a panel subjected to inplane loading. 


For this configuration the geometry parameter is given as 


a ni HL 
b V 1 + c2 


(34) 


and the corresponding expression for the flutter parameter is 


FP = 


D, f(N) 


1 




(35) 


24 




FIGURE 5 PANEL WITH INPLANE LOADING 


The quantity equations (34) and (35), is the ratio of the inplane load 
N-j in the x direction to the critical buckling load Thus 


Pcp - (35) 

As indicated in Section 3.0 by equation (13), the basic form for 
GP including the influence of inplane loads is based on the two mode 
flutter solution. The C /K^ term in equation (34) was added to yield 
results consistent with equation (20). The form of the flutter parameter 
given by equation (35) was obtained by matching trends in experimental panel 
flutter data, primarily from reference 12. 

To make flutter predictions for a panel with inplane loading, the magni- 
tude of is required for use in equation (36). Evaluation of this buckling 
load parameter must account for both the and loadings (fig. 5). Procedures 
for obtaining are presented in references 13 and 14. 
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Often in preliminary design neither the magnitude of the inplane loads 
nor the critical buckling load are known with any degree of confidence. As 
indicated in reference 15, for a given panel configuration subjected to in- 
plane loading, the minimum dynamic pressure for panel flutter occurs at the 
transition point between an unbuckled and buckled panel. Thus, a conservative 
approach, when the effects of inplane loading are desired, would be to assume 
a value for of one. Equation (34) shows that when equals one the 
geometry parameter equals zero. In this situation the magnitude of the flutter 
parameter may be obtained from equation (14). The flutter critical dynamic 
pressure is then obtained from equation (35). 

The design approach is limited to panel configurations with normal in- 
plane loads and does not account for the possibility of inplane shear loading. 

A recent theoretical investigation of this question of combined normal and 
shear inplane loads is presented in reference 16. The conclusion of this study 
is that designs for panels at buckling due to inplane shear or combinations of 
inplane shear and normal loads will be conservative if it is assumed that the 
panels buckle due to N-j alone. 

With the established value of FP, the magnitude of q/f(M) is determined 
from equation (35). The flutter critical dynamic pressure q is then obtained 
after consideration of the effects of static pressure differential. Section 
4.4; structural damping. Section 4.5; and Mach number. Section 4.6. 

4.4 Static Pressure Differential - It has been shown experimentally 
(refs. 17 and 18) that a static pressure differential Ap across the panel 
has a significant effect on the panel flutter boundary. This effect is due 
to the inplane stresses that are induced by this pressure differential. 

These stresses are always tensile regardless of the direction in which Ap 
acts and increase the effective stiffness of the panel. Thus, the effect of 
this differential pressure is to raise the panel's flutter boundary. 

In reference 17, the effect of a pressure differential is described 
through a nondimensional parameter given as 

„ Ap a^ (37) 

^ D h 
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where an isotropic panel is assumed. For an orthotropic panel an analogous 
nondimensional parameter of the form 


P 


f 


Ap a 

OXn 

1 eq 


(38) 


has been defined. In this expression h^^ is an equivalent isotropic panel 
thickness quantity related to the orthotropic panel stiffness properties 
and is defined as 


h 


_ 1 
eq “ 2“ 



+ 



(39) 


The dynamic pressure parameter including the influence of a static 
pressure differential becomes 



where Qp is the ratio of the dynamic pressure parameter qp/f(M), which 
includes the effects of Ap, to q/f(M) which does not include pressure 
differential effects. Shown in figure 6, as a function of P^, is the 
parameter Qp for panels with several length to width ratios. These trends 
are based on experimental data discussed in references 17 and 18. 

The results summarized in figure 6 form the basis for including the 
influence of a Ap in the panel design. Following the determination of the 
flutter critical parameter q/f(M), as covered in the preceding sections, 
the magnitude of is calculated from equation (38). With this value of 
P^, the appropriate Qp is obtained from figure 6 and the updated quantity 
qp/f(M) obtained. 

4.5 Structural Damping - The influence of structural damping in a 
panel system with inplane loading is accounted for through a damping factor 
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D^, defined as 


(41) 


[ 


1 + 100 P 


CR ^9 


In the above is the buckling load ratio of equation (36) and the quantity 
Ag is expressed as 


Ag = g - 0.01 


(42) 


where g is the damping coefficient assumed in the system. 


With the definition of D^, equation (41), the dynamic pressure parameter 
including the influence of damping is obtained from the following relationship. 

The quantity q/f(M) in equation (43) is the result of the flutter prediction 
employing the flutter-free design boundary as discussed in the previous sec- 
tions, and q^/f(M) is the updated quantity including the effects of damping. 

The quantity Ag is defined as indicated by equation (42) because appli- 
cation of the design procedures produce results which compare well with 
theoretical results for panels with inplane loading and panel damping repre- 
sented by a structural damping coefficient of 0.01. Thus, results obtained 
from the design boundary may be thought of as including "nominal" damping 
of this magnitude. The parameter Ag is introduced in the damping factor to 
account for panel damping which is higher than this nominal value. 

The above form of 0^ was obtained through comparison with the results 
of references 15 and 19. These references present trends, obtained during 
detailed theoretical investigations, of flutter critical flow conditions as 
influenced by inplane loading and structural damping. 
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4.6 Mach Number Effect - Use of the dynamic pressure parameter q/f(M) 
to describe the flutter critical free stream flow conditions is based on the 
commonly used linear piston theory, reference 20, to represent the aerodynamic 
forces acting on the panel. As indicated in reference 20, experience has 
shown that use of this aerodynamic theory in conjunction with f(M) defined 
as 


f(M) =-\f M^-1 


(44) 


yields flutter boundaries that exhibit excellent agreement with those given 
by more refined theories, provided the Mach number is higher than approxi- 
mately 1.6 to 2.0. 

For application with the flutter- free panel design boundary, f(M) is 
assumed to take the form of equation (44) for Mach numbers greater than 2.0 
and is defined as shown in figure 7 for Mach numbers in the range of 1.0 to 
2.0. These curves were derived from experimental data discussed in refer- 
ence 18. Following determination of the parameter q/f(M) through application 
of the design boundary and including modifications for potential static pres- 
sure differential and damping effects, the flutter critical free stream 
dynamic pressure is obtained for the Mach number correction factor defined by 
either equation (44) or figure 7. 
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5.0 EXAMPLE APPLICATIONS OF DESIGN PROCEDURE 


Several examples are presented which illustrate the use of the panel 
design procedures detailed in Section 4.0. Both isotropic and orthotropic 
panels with various combinations of support conditions, inplane loading, and 
flow angularity are discussed. In all cases, it has been assumed that the 
applicable factors of safety have been considered in arriving at the stated 
free stream dynamic pressure requirements. 

5.1 U niform Panel with Inplane Loads - It is desired to determine the 
thickness of a uniform aluminum panel required to preclude flutter when the 
panel is subjected to an inplane load that is 75 percent of its buckling load. 
A sketch of the panel is shown in figure 8 and details of this configuration 
are as follows: 

Simply supported panel 
a = 0.25 m 

b = 0.65 m 

E =69 GN/m^ 

PCR =0.75 

M = 2.5 

q =40 kPa 

Figure 4 shows that the design approach for a panel with inplane 
loading is detailed in Section 4.3. For an isotropic panel, the panel stiff- 
ness properties D.j and D .|2 are equal and equation (34) becomes 

The above expression also reflects that for a panel simply supported at all 
edges, the quantity approaches infinity and thus, the C /K^ term of 
equation (34) approaches zero. Thus the magnitude of the geometry parameter 
for this particular panel configuration becomes: 

GP = 0.192 
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FIGURE 8 UNIFORM PANEL WITH INPLANE LOADING 


From equation (14) with an of 0.027 or from figure 3, the magnitude 
of the flutter parameter is determined to be 

FP = 0.00532 

From equation (35) the panel stiffness parameter to preclude flutter may be 
expressed as 

where q is the free stream dynamic pressure requirement. Note that it has 
been assumed that no modifications to the dynamic pressure which account 
for damping or pressure differential effects are required. The Mach number 
correction factor (Sec. 4.6) is given as 

f(M) =\/m^-1 = 2.29 

and the resulting panel stiffness requirement becomes 

D = 88.26 N-m 
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For an isotropic panel, the panel stiffness is expressed as 


where E is the material modulus of elasticity, h the panel thickness, and v 
the material Poisson's ratio. Using a value of 0.3 for v, the panel thick- 
ness requirement to preclude flutter becomes 

h = 2.4 mm 

5.2 Elastically Supported Orthotropic Panel - ftn orthotropic panel 
simply supported at the leading and trailing edges and elastically supported 
at the stream edges as shown in figure 9 is to be evaluated for flutter. 

The influence of stream edge flexibility on the flutter critical dynamic 
pressure at a Mach number of 1.7 is required to aid in the definition of the 
support springs. The physical characteristics for the panel are as follows: 


a = 0.45 m 

b = 0.90 m 

E =70 GN/m^ 
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= 10 N-m 
Dg = 2000 N-m 
0^2 = 750 N-m 
Ap =2 kPa 

Kp = range of 2 to 20 MN/m/m 
Zero flow angle 


From figure 4 the design procedure to be followed in this instance is 
detailed in Section 4.1. The expressions for the geometry and flutter 
parameters are given by equations (20) and (21): 



'^ 12/°1 
1 + C^/Kj 


FP = 


0^ f(M) 
q a^ 


The fact that the leading and trailing edges have complete deflectional 
restraint (simply supported) has been reflected in the above by setting the 
nondimensional support parameter [eq. (22)] equal to one. Thus, for 
the particular panel configuration of interest 

75 

1 + 28.I3/K3 

q 109.7 
flMT FP 



It is assumed that the stream edge supports take the form of a running 
spring. From figure B-2 of Appendix B the nondimensional support stiffness 
is given as 


K 


S 


3 

TT 


b 


D 


3 


2 
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which becomes 


K3 = 1 .18 X 10'^ Kp 

for the configuration under consideration. 

From figure 7, with an a over b ratio of one half and a Mach number 
of 1.7, the Mach number correction factor has a magnitude of 


f(M) = 1.23 


From equation (38), the pressure differential parameter is 


P 


f 


Ap a 


4 


eq 


where h^^ is given by equation (39) and has the following value for this 
particular configuration 


h = 4.1 mm 
eq 

Thus, the pressure differential parameter becomes 

Pf = 2000 

and from figure 6 the pressure differential correction factor is determined 
to be 

Qp = 1.2 

Combining the above, the relationship between the flutter critical dynamic 
pressure and the flutter parameter is given as 

q — Qp = __ — 
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The resulting influence of the stream edge support stiffness on the 
critical dynamic pressure is illustrated in figure 10. It must be noted 
that since GP is always less than five, the corresponding values for FP are 
obtained from equation (14). The zero moment boundary condition value 
of is employed for the calculations summarized in figure 10. 


5.3 Orthotropic Panel with Flow Angularity - The flutter sensitivity 
of a clamped orthotropic panel to angular flow is to be evaluated. Physical 
characteristics of the panel in question, shown in figure 11, are as follows: 
a = 1 .1 m 

b = 0.7 m 

= 2500 N-m 
= 10 N-m 

d-^2 “ 

M = 3.5 


From figure 4 the procedures of interest for this case are presented 
in Section 4.2. For a panel with complete deflectional restraint (fully 
clamped) the desired relationships for GP and FP are given by equations (32) 
and (33). These relationships become 


and 


:yfc 


GP = 0.222\/Cos^ A + 41 Sin^ A 


q = 


6300 


FP [Cos^^ A + 64.4 Sin^ a] 


for this particular panel configuration. Reflected in the above expression 
for the flutter critical q is the evaluation of the Mach number correction 
factor from equation (44). 

f(M) =Vm^-1 = 3.354 
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FLUTTER CRITICAL DYNAMIC PRESSURE - kPa 













The preceding expression for GP has been evaluated for various magnitudes 
of flov/ angularity A. The corresponding value for FP was obtained from 
equation (14) using the zero slope boundary condition value for E^. The flutter 
critical dynamic pressure was then obtained and the resulting relationship 
between q and A is shown in figure 12. 



q 


FIGURE 11 ORTHOTROPIC PANEL WITH FLOW ANGULARITY 
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FLUHER CRITICAL DYNAMIC PRESSURE - kPa 
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FIGURE 12 FLUTTER CRITICAL DYNAMIC PRESSURE VERSUS 
FLOW ANGULARITY 
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6.0 DESIGN BOUNDARY VERIFICATION 


The preceding section describes a design approach, in conjunction with 
the flutter- free panel design boundary, developed for use in the preliminary 
design of flutter free lightweight external panels. In this section, the 
validity of the approach, as it applies to various panel configurations, 
is evaluated by comparing results from the preliminary design predictions 
with appropriate existing experimental data. Where no experimental data 
exist, comparisons are made with the results of detailed theoretical analyses. 
The experimental panel flutter data employed in the following sections are 
tabulated in Appendix D and cross-referenced to the appropriate figures in 
the following sections. 

6.1 Basic Panel - A large amount of experimental and flight flutter 
data for the basic panel configuration is available. In this context, "basic" 
panel configuration refers to a flat orthotropic panel (fig. 1) with either 
complete edge deflectional restraint or edge support flexibility effects. 

The influence of additional parameters such as flow angularity, inplane loads, 
etc., is considered in later sections. 

The panel flutter data for the basic panel configuration obtained from 
references 10, 18, and 21 are shown in figure 13. The data of references 18 
and 21 are for isotropic panels with total edge deflection restraints. As 
seen in figure 13, a majority of these data compare well with the flutter-free 
panel design boundary. The reference 18 data were obtained for a number of uni- 
form panels of the same thickness and length to width ratios. For this panel 
configuration, the 37 experimental data points (represented by the straight 
line connecting the circular symbols) all fall within the flutter boundary. 

The data of reference 21 are for clamped isotropic panels of varying thickness 
and length to width ratios. These experimental results also compare well with 
the design boundary. 

Data from an experimental investigation on the flutter characteristics 
of orthotropic panels are given in reference 10. These results are for a 
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FIGURE 13 FLUTTER DATA FOR BASIC PANEL CONFIGURATIONS 


wide variety of panel configurations and stream edge flexibility support 
conditions. Their comparison with the flutter-free panel design boundary 
shows the boundary to be conservative. Indicated in figure 13 is the importance 
of the C /Kg term in the geometry parameter of equation (20). The solid points 
shown in this figure present the reference 10 data without accounting for the 
stream edge support flexibility while the open symbols include the influence 
of this support flexibility. It can be seen that when the edge flexibility 
effects are neglected the design boundary is very unconservative. 

The lack of experimental data for panels with leading and trailing edge 
support flexibility precludes the comparison of the design boundary with 
test data. As an alternate, predictions obtained through application of the 
flutter- free panel design boundary are compared with theoretical results of 
reference 9. This comparison is shown in figure 14 and as can be seen there 
is adequate correlation between the trends obtained with the design boundary 
and the theoretical results of reference 9. 

6.2 Flow Angularity - A limited amount of experimental flutter data, 
as presented in references 8, 11, and 22, exist for panels exposed to flow 
not parallel to one 'side. Note that the results presented in reference 22 
are a duplication of the reference 11 information. All of these data, except 
for panel 1 of reference 11, are for orthotropic panels with leading and trail- 
ing edge support flexibility with respect to zero angle flow. The panel 1 
configuration was an orthotropic panel clamped along all sides. 

Comparison of these data with the design boundary is shown in figure 
15. As can be seen, for a majority of these experimental data the correlation 
with the design boundary is conservative. Shown in figure 16 are the experi- 
mental and predicted x's as a function of flow angle for the panel 2 (ref. 11) 
data. This figure illustrates the adequacy of the design approach in pre- 
dicting the influence of flow angle on the flutter critical flow conditions 
of an orthotropic panel. 
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U/ Dj 

FIGURE 14 CORRELATION WITH REFERENCE 9 THEORETICAL RESULTS 


6.3 Inplane Loads - Experimental data obtained during an extensive 
wind tunnel test program conducted to evaluate the influence of inplane loads 
on panel flutter (ref, 12) are compared with the design boundary in figure 17. 
These data are for isotropic panels with a wide range of panel a over b ratios. 
During the test program mechanical inplane loading was introduced by hydraulic 
actuators. As indicated in this figure, the correlation between these data 
and the flutter-free panel design boundary is conservative. 
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FIGURE 15 PANEL FLUTTER DATA INCLUDING INFLUENCE OF FLOW ANGULARITY 
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FIGURE 16 CORRELATION WITH REFERENCE 11 EXPERIMENTAL RESULTS 
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FLUTTER PARAMETER FP 

FIGURE 17 FLUTTER DATA FOR PANELS WITH VARIOUS a OVER b RATIOS AND INPLANE LOADS 


It is assumed in reference 12 that the test panels have complete 
deflectional support at the edges. However a review of the test setup indi- 
cates that some stream edge support flexibility existed. This effect, as 
presented in Table D-6 was included in the calculation of GP for the data 
shown in figure 17. 

Additional experimental panel flutter results, including the influence 
of inplane loads, are given in references 7, 19, and 23 through 26. Test data 
as presented in these references are plotted on the panel design boundary in 
figure 18 The overall comparison between these data and the design boundary 
is in general conservative. For all these panels the inplane loads are the 
result of panel aerodynamic heating. The data shown on figure 18 are for 
isotropic panels except those discussed in references 7 and 26. 

6.4 Static Pressure Differential - The design procedure to be followed 
in evaluating the influence on panel flutter of a pressure differential is 
based on experimental results presented in references 17 and 18. These wind 
tunnel data were obtained for isotropic panels having complete restraint 
(clamped or simply supported) at the panel boundaries. 

Limited amounts of additional data including pressure differential ef- 
fects are available for correlation with the recommended design approach. 

The available data (ref. 10) are compared with the design boundary in figure 
19. These data are for a highly orthotropic panel with stream edge support 
flexibility and an a over b ratio of one. 

The closed symbols are for this experimental data before the pressure 
differential correction factor Qp has been applied. Application of the 
procedures of Section 4.4 to modify the experimental q/f(M) value by the 
appropriate magnitude of Qp results in the comparison with the design bound- 
ary as indicated by the open symbols in figure 19. 
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FIGURE 18 PANEL FLUTTER DATA INCLUDING INFLUENCE OF INPLANE LOADS 




GEOMETRY PARAMETER 



10 '® 10 '^ 10 '^ 10 '^ 10 '^ 

FLUTTER PARAMETER - FP 

FIGURE 19 PANEL FLUTTER DATA INCLUDING THE INFLUENCE OF STATIC PRESSURE DIFFERENTIAL 



other data with pressure differential effects included are available 
in reference 10. However, these configurations result in magnitudes for the 
pressure differential parameter which exceed the range covered in figure 6. 
This indicates that the design procedure of Section 4.4 is limited to panel 
configurations whose characteristics result in magnitudes of covered by 
figure 6. 


6.5 Structural Damping - A procedure is described in Section 4.5 
which allows eval uation of the influence of structural damping in combination 
with inplane loads. With this procedure the dynamic pressure parameter q/f(M) 
obtained through application of the design boundary is modified by the damping 
factor Df [eq. (41)] to account for the presence of the damping. 

For a panel with inplane loading, results obtained from application of 
the flutter- free panel design boundary with no added damping compare well with 
detailed theoretical analyses that include a structural damping coefficient 
of 0.01. Such a comparison is shown in figure 20, where both experimental 
data and theoretical results for a clamped panel detailed in reference 19 are 
shown. Thus, results obtained from the design boundary may be thought of as 
including "nominal" structural damping on the order of 0.01. For panels 
having higher damping characteristics (g > 0.01) and inplane loads, the damp- 
ing factor 0^ as defined by equation (41) is to be employed. 

Many data points for panels with inplane loading shown in figure 17 
appear quite conservative. As discussed in reference 12, significant struct- 
ural damping, up to five percent, was measured for a number of the panels. 

The presence of such magnitudes of damping can have a large influence on the 
flutter susceptibility of panels with inplane loading. Referring to equation 
(41), the damping factor is dependent on the magnitudes of both the structural 
damping and inplane loading. Figure 21 shows the influence of structural 
damping on panel 4 of reference 12. The lower flutter boundary results from 
the design boundary (g = 0.01) and the upper boundary results from the damping 
factor defined by equation (41). The good agreement between the experi- 
mental data and the predicted boundary for the average measured structural 
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FIGURE 20 CORRELATION WITH REFERENCE 19 CLAMPED PANEL RESULTS 
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FIGURE 21 INFLUENCE OF HIGH STRUCTURAL DAMPING FOR A PANEL 
WITH INPLANE LOADING 





damping g = 0.028 indicates that equation (41) is sufficiently accurate for 
flutter design of panels with significant structural damping. 

6.6 Mach Number Effect - The definition of the Mach number correction 
factor f(M) as given in figure 7 is based on experimental panel flutter 
results in the low supersonic Mach number regime. The trends illustrated in 
figure 7 are established in reference 18 and are based on results obtained 
during wind tunnel test programs such as those discussed in references 18 
and 21. 


Further verification of these relationships between f(M) and Mach 
number is illustrated in figure 22. Here data obtained during extensive 
panel flutter test programs detailed in references 27 and 28 are compared 
with predictions obtained with the flutter- free panel design boundary when 
modified by f(M). The data from reference 27 are for a clamped isotropic 
panel while those of reference 28 were obtained for a panel having less than 
a fully clamped boundary condition. The spread in test data indicated in 
figure 22 for a particular Mach number is related to the boundary layer 
thickness at the test conditions. 
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7.0 CONCLUDING REMARKS 




Criteria and procedures for their application in the preliminary design 
of flutter free lightweight external surface panels have been presented in a 
format which permits their use by individuals without prior panel flutter ex- 
perience. The major area of application for these procedures is in the design 
of orthotropic metallic panels which stand off from primary structure. Examples 
which illustrate the use of these criteria and procedures are given. 

The panel design approach presented in this report is a logical extension 
of previous work done in this area. The design criteria account for most of 
the interacting parameters which significantly influence panel susceptibility 
to flutter. Parameters included in these design procedures are 
0 Panel aspect ratio 
0 Panel orthotropic properties 
0 Support flexibility 
0 Inplane loads 
0 Static pressure differential 
0 Flow angularity 
0 Structural damping 

The validity of the flutter-free panel design boundary (fig. 3) and its 
application (Section 4.0) to preliminary panel design has been demonstrated 
by comparison of predicted flutter results with experimental data. Included 
in this verification process are data for panels ranging from isotropic to 
those having highly orthotropic stiffness characteristics. In addition, data 
for panels including the influence of parameters such as edge support condi- 
tions, inplane loads, flow angularity, etc., were available and used. As dis- 
cussed in Section 6.0, application of the design procedures resulted in the 
conservative prediction of flutter points for a vast majority of the existing 
experimental data. 

External panels on discrete flexible supports are being strongly considered 
for use as a high performance vehicle thermal protection system. The influence 
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of this type of support condition was not incorporated into the design proce- 
dure due to a limited availability of information (ref. 29) for this configura- 
tion. Further theoretical considerations and corresponding panel flutter test 
programs are needed to investigate this support condition. 

Limited experimental data exist for the flutter of orthotropic panels 
including the influence of flow angularity. In addition, no data exist for 
panels subjected to the combined conditions of inplane loads and flow angular- 
ity. Comprehensive theoretical investigations of these aspects of the panel 
flutter problem have been conducted and are documented in the literature. 
Further testing to provide the data needed to evaluate these theoretical 
results is required. 
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APPENDIX A 

ORTHOTROPIC PANEL STIFFNESS PARAMETERS 

Of basic importance to the use of the flutter-free panel design boundary 
is knowledge of the panel orthotropic stiffness properties. A number of 
techniques which may be employed to obtain these quantities are presented in 
this appendix. 


For panels with closed corrugations [fig. A-l(a)] the use of the proce- 
dures detailed in reference 30 are recommended. Determination of the panel 
stiffness properties for the other two configurations shown in figure A-1 is 
discussed in reference 31. For a panel with equidistant ribs [fig. A-l(b)] 
reference 31 gives the following expressions for D.j , D 2 , and D -|2 

0 = i 

’ 12(l-v^) [Jl-t + (h/H)^t] 


D 


2 



D 


12 ■ 


Eh^ , GJ 
12(l-v^) 


(A-1) 


where I^ is the moment of inertia of the cross-hatched area defined in figure 
A-l(b) and GJ is the torsional rigidity of one rib. 


The stiffness properties for the corrugated sheet [fig. A-l(c)] as de- 
fined in reference 31 are expressed as 

^ ^ 12(1-v2) 

(A- 

D2 = Ell 


'12 


_ s Eh’’ 

1 12 (l+v) 
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(a) CLOSED CORRUGATIONS 


T f-iH 


(b) EQUIDISTANT RIBS 


z 



(c) CORRUGATED SHEET 


FIGURE A-1 ORTHOTROPIC PANEL CONFIGURATIONS 


where 




(A-3) 


and 


= 

1 2 


1 


.81 ' 

1 + .625 (H/Jl)^ . 


(A-4) 


An approach for obtaining initial estimates for the stiffness properties 
of an orthotropic panel is implied by Timoshenko in reference 31. With this 
approach, the panel bending rigidities are related to the properties of a beam 
of unit width. Thus, considering the case illustrated in figure A-1(b), we 
have 


D 



0 


2 



(A-5) 


where I.| has previously been defined. Following the approach mentioned in 
reference 31, the quantity 0^2 estimated as 

“l2 ■ 

The above is not a recommended procedure for general usage, but does allow 
initial estimates of the panel's stiffness characteristics to be readily made. 


Experimental procedures for obtaining these panel parameters are discussed 
in references 32 and 33. These experimental approaches do not have application 
during preliminary design, but are mentioned here for completeness. 
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APPENDIX B 

EDGE SUPPORT STIFFNESS EXPRESSIONS 

Techniques have been developed for the estimation of the nondimensional 
stiffness parameters for various panel edge support configurations. An energy 
approach is used to obtain these expressions. In this approach, for an assumed 
deflection pattern, the strain energy stored in an idealized support spring 
configuration (fig. B-1) is equated to the energy stored in the configuration 
of interest. 

A summary of these results is given in figure B-2. Expressions are pre- 
sented in this figure for the support stiffnesses along with sketches defining 
each support's geometric characteristics. The results presented in figure B-2 
are for support springs located along the panel stream edges as illustrated in 
figure B-1. Similar expressions for leading and trailing edge support springs 
are obtained by interchanging the geometric quantities associated with the 
panel coordinate system. The stiffness expressions are given in the nondimen- 
sional form required for application with' the flutter-free panel design bound- 
ary. 



FIGURE B-1 IDEALIZED STREAM EDGE SPRING SUPPORTS 


65 



SUPPORT 

CONDITION 

SUPPORT 

TYPE 

SUPPORT CONFIGURATION 

NONDIMENSIONAL SUPPORT STIFFNESS 

1 

RUNNING 

SPRING 


h 

h = ^ 

■n^ D 2 

2 

DISCRETE 

SPRING 



3 

EDGE 

SUPPORT 

PANEL 



4 

TAPERED 
EDGE SUPPORT 
PANEL 


WHERE 

5 

DOUBLER 

Vei^ 

(Ifd) 


FIGURE B-2 NONDIMENSIONAL STREAM EDGE SUPPORT STIFFNESS EXPRESSIONS 
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A summary of the stiffness properties for a number of support clips or 
standoffs is presented in figure B-3. These quantities are presented in terms 
of the stiffness associated with a discrete spring. The corresponding running 
or line spring for a number of discrete springs is defined as 


K 


D 



(B-1) 


where N is the number of discrete springs located down the panel side of 
length a. 


To illustrate the procedure that was employed to obtain the information 
presented in figure B-2, the analysis of an edge support panel (support condi- 
tion 3) will be presented in detail. Referring to figure B-1, the strain 
energy stored in the idealized support springs is given as 




S"' 


dx 



(B-2) 


where it has been assumed that the deflection at the support springs is given 
as 


w = C sin — 


Introducing the- nondimensional spring rate notation 

3 


Kn = 




equation (B-2) becomes 


D 3 n 

7T 


u, ■ c2 f Dj K„ 


(B-3) 


(B-4) 


(B-5) 


Since for stream edge support flexibility we have 


K3 = K, 


(B-6) 
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MAIN PANEL 


MAIN PANEL 



FIGURE B-3 STIFFNESSES OF PANEL EDGE SUPPORT SPRINGS 





equation (B-5) may be written as 



2 a \ 3 
4 


Dg Kg 


(B-7) 


For the edge support panel illustrated in figure B-4 the energy stored in 
the edge panel is given by 











(B-8) 


It is assumed that the deflection throughout the edge panel is given as 




Sin 


TtX 

a 


(B-9) 


The use of this deflection shape is analogous to assuming that the main panel 
is infinitely stiff relative to the stiffness of the edge panel in the y direc- 
tion. Combining equations (B-8) and (B-9) yields the following expression for 
the strain energy stored in the edge panel 


1 r2 


= -L c'- n r 82. rii** 4. / If X T5 I 1 

4 ^ ^0 4 I- ^ ^ ^ ( il ^ + 12 3 


68 /tt\4 48 


(B-10) 


Equating equations (B-10) and (B-7), the expression for the nondimensional 
stiffness of an edge support panel becomes 


K 


S 





17 

35 


(f) 


12 (^)5 
5 'aai 



(B-11) 
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APPENDIX C 

PANEL ON ELASTIC FOUNDATION 


The question of the flutter sensitivity of a panel on an elastic founda- 
tion falls into two categories. The first of these, with a rigid substructure, 
is illustrated in figure C-l(a). It consists of a flexible panel exposed to 
the flow while supported by an elastic medium or foundation. This elastic 
foundation is in turn mounted to an essentially rigid base. In the second, 
the lower surface of the elastic medium is mounted to a second flexible panel 
or elastic substructure as shown in figure C-l(b). These two different con- 
figurations are treated separately in the following sections. 

C.l Rigid Substructure - In the absence of damping it can be shown 
(ref. 34) that the critical flow conditions for flutter of a simply supported 
isotropic panel resting on an elastic foundation with a rigid substructure 
are unaffected by the foundation. In reference 34 a Winkler foundation model 
was employed wherein a point load on the foundation gives rise to a deflection 
at the same point but nowhere else. The distinction between this Winkler 
foundation representation and a foundation model which allows deflections at 
points surrounding a point load is discussed in reference 35. 

The fact that a Winkler foundation has no influence on the flutter of a 
simply supported panel, can be explained by the way in which the foundation 
contributes to the panel dynamic characteristics. The effect of the founda- 
tion is to increase the panel natural frequencies, but not to change the fre- 
quency separation between the various modes. This results in the same flutter 
q for a given panel configuration with or without the elastic foundation. As 
discussed in reference 34, when aerodynamic damping or an effective structural 
damping (viscous type) is included in the analysis the flutter conditions can 
be significantly affected by the foundation depending upon the magnitude of 
the foundation stiffness. 
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(a) Rigid Substructure 



FIGURE C-1 PANEL ON ELASTIC FOUNDATION 



Results of a two mode flutter solution for an isotropic panel presented 
in reference 34 have been extended to the case of an orthotropic panel rest- 


elastic 

foundation. 

. These 

results 

are expressed as 

A = 

18.26 

15 + 6 

^^12 

Dl 



I 1 

+ 

O 

^12 / 
°1 \ 

f/ - 



where A is defined by equation (6). In this expression, g.^- is a damping co- 
efficient which combines both aerodynamic and structural damping, K is the 
elastic foundation stiffness parameter given as 





with other quantities defined in figure C-l(a). 

The design procedure for this panel configuration employs the flutter- 
free panel design boundary to estimate a flutter critical condition neglecting 
gy and elastic foundation effects. With equation (C-1), the ratio of flow 
parameters with and without elastic foundation and damping effects is obtained. 
The critical parameter obtained from the design boundary is then modified by 
this ratio to account for damping and foundation stiffness. 

Results obtained from this approach for a panel configuration with a rigid 
substructure are shown in figure C-2. The relationship between the flow para- 
meter ratio A/A^ obtained from equation (C-1) as a function of the nondimen- 
sional foundation parameter K [eq. (C-2)] and damping is illustrated. In this 
context Aq is the flow parameter at flutter neglecting damping and elastic 
foundation effects while A is the corresponding parameter including these ef- 
fects. These results are compared with the analytical results from reference 
34. 
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FLOW PARAMETER RATIO 



FIGURE C-2 COMPARISON OF TWO MODE SOLUTION INCLUDING ELASTIC 
FOUNDATION EFFECTS WITH REFERENCE 34 RESULTS 




As indicated by the data presented in figure C-2, this procedure has ap- 
plication only at lower values of K (K _< 1000). It must be emphasized that 
this approach has the inherent limitations of a two mode flutter solution and 
a reliable panel design effort must include a detailed flutter analysis includ- 
ing the influence of the elastic foundation. 

C.2 Elastic Substructure - Theoretical considerations for two parallel 
isotropic simply supported elastic panels interconnected by an elastic medium 
are given in references 36 and 37. Results shown in these references illus 
trate the complex interplay between the stiffness characteristics of the elas- 
tic medium, panel configurations, and magnitude of inplane loading. Represen- 
tation of these complex interactions in terms of design criteria does not 
appear feasible. However, the two mode theoretical approach of references 36 
and 37 has been modified for orthotropic panels and additional panel flutter 
trends obtained. 

Following the analytical procedure of references 36 and 37 and neglecting 
the influence of inplane loads, the two mode solution for parallel, elastically 
interconnected simply supported orthotropic panels is 
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where ^ is the mass per unit area for the particular panel. In figure C-l(b) 
the superscripts U and L refer to the upper and lower panels, respectively. 
The form of equation (C-3) assumes deflection solutions for the upper panel 
as follows 


(x.y.t) = wl/jt) sin f- sin W^^(t) sin 


2ttX 

a 




(C-4) 


with similar expressions for the lower panel 


Equation (C-3) is nondimensional ized with respect to the lower panel by 
employing the following relationships 


0l ■ ' 


-U _ -L 
p “ P 


Thus, equation (C-3) may be written as 


X = 


K = 


2qa' 


f(M)D 


k a^ 
4n L 

77 


(C-5) 


TT D.J 

-LA 
p a 


[1 - 2(^) ' . K 


8 X 

I -T 


-K 

0 


8 

T "T 


-K 


[16 + 8 (f) 2 + (Ji) ^ Q 

0 [1 + 2 {^) '■ (f) ^ + (^)'- (f)'^ ] + K 


-K 


[16 + 8 + (^)4f)^] + K 


W 

(C-6) 


n 0 0 

0 n 0 

0 0 1 

_ 0 0 0 


0 

0 

0 

1 , 


W = 


0 
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The eigenvalue problem associated with equation (C-6) was formulated and 
programmed on the digital computer. An interative -procedure was then used to 
obtain the flutter critical x for a given panel configuration. In this pro- 
cedure the value of x is systematically varied and the lowest value of x for 
which two of the eigenvalues coalesce is sought. The procedure was used to 
generate trends for flutter of coupled orthotropic panels. These trends are 
presented so that the designer is aware of the potential flutter problems with 
parallel elastically coupled panels. 

Typical results for interconnected parallel panels are illustrated in 
figure C-3 where the critical flow parameter x as a function of panel a over b 
ratio and elastic medium stiffness is shown for identical upper and lower iso- 
tropic panels (e = n = 1.0). As can be seen there is a significant dip in 
the critical flow parameter as a function of the a over b ratio. A plot of 
the elastic medium stiffness versus panel aspect ratio for this dip in critical 
X is given in figure C-4. 


A portion of the results presented in figure C-3 can be compared with 
the data presented in reference 37. Shown on figure C-3 is the comparison 
which indicates that the trends presented here are consistent with the refer- 
ence 37 results. 


If it is assumed that the two panels have different thicknesses, but are 
constructed of the same material, the quantities e and n of equation (C-5) 
become 



n = T 


r 


(C-7) 


where T^ is the thickness ratio expressed as 

T^ = h^/h*- 


(C-8) 
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PANEL LENGTH TO WIDTH RATIO - a/b 


FIGURE C-4 FLUTTER CRITICAL FOUNDATION STIFFNESS PARAMETER AS A 

FUNCTION OF LENGTH TO WIDTH RATIO FOR EQUAL ISOTROPIC PANELS 


Two mode results for an upper panel having twice and half the thickness of the 
lower panel are compared with results for equal panel thicknesses in figure c-5. 
These configurations have a critical foundation stiffness parameter K values 
which result in low values for the flow parameter at flutter. 

An analysis of elastically interconnected panels representative of the 
space shuttle thermal protection system is presented in reference 38. In this 
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FLOW PARAMETER - X/ir 



FIGURE C-5 FLUTTER RESULTS FOR PARALLEL ELASTICALLY INTERCONNECTED 
ISOTROPIC PANELS HAVING VARYING THICKNESS RATIOS 






configuration the upper panel consists of a relatively thick ceramic tile. 
Included in this investigation is the influence of several individual tiles 
supported by a common lower panel and interaction between adjacent tiles. The 
designer is referred to this work for detailed discussions of this type of 
panel configuration. 

In the absence of more extensive theoretical and/or experimental results, 
the above two mods approach allows the evaluation of trends for parallel panel 
configurations. These procedures could be employed during preliminary design, 
but as with the case of a rigid substructure, detailed flutter studies must be 
conducted during later design efforts. 
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APPENDIX D 

PANEL FLUTTER TEST DATA 


A tabulation of the experimental panel flutter data used in the verifi- 
cation of the flutter-free panel design boundary. Section 6.0, is presented 
on the following pages. In each case, the panel configuration (isotropic or 
orthotropic, flow angularity, inplane loads, etc.) is indicated. A separate 
table is given for the data obtained from a particular reference for a given 
panel configuration. Also presented is a cross reference with the figure 
number of Section 6.0 in which each set of data is plotted. 
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TABLE D-1 


PANEL FLUTTER DATA FROM REFERENCE 10 
Orthotropic panel with steam edge support flexibility 
Data shown in figure 13 


Panel 

No. 

Support 

Condition 


a 

(m) 

U-1 

3 

5.0 

0.61 

U-2 

3 

2.5 

.61 

V-1 

3 

3.7 

.61 

V-4 

3 

3.4 

.61 

H-1 

3 

.8 

.61 

S-1 

4 

1 .5 

.946 

S-2 

4 

.7 

.946 

S-3 

4 

.7 

.946 

S-4 

4 

.7 

.946 

S-5 

4 

1 .5 

.946 

S-6 

4 

1 .5 

.946 

S-7 

4 

1 .5 

.946 

S-8 

4 

.7 

.946 

S-9 

4 

.7 

.946 

S-IO 

4 

.7 

.946 


0 

(m) 

•Jl 

(N-m) 

U2 

(N-m) 

(N-m) 

0.61 

0.35 

204 

52 

.61 

1 .03 

1640 

376 

.61 

1.9 

2180 

344 

.61 

5.4 

5420 

1190 

.61 

.94 

4150 

798 

.441 

1 .15 

2700 

856 

.408 

1 .19 

5350 

1760 

.408 

1.25 

5060 

1800 

.408 

1 .33 

4790 

1840 

.441 

1 .02 

3060 

800 

.441 

.33 

987 

258 

.441 

.37 

870 

277 

.408 

-.94 

6800 

1520 

.4C8 

.31 

2200 

490 

.408 

.47 

5460 

1210 


Note: Support condition definition and Kg calculation based on 

figure B-2, Appendix B and 


X 


_ 2q.a3 
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TABLE D-1 (Concluded) 


PANEL FLUTTER DATA FROM REFERENCE 10 


Panel 

M 

q 

Qp(l) 

Gp(2) 


No. 


(k Pa) 



(xlO-5) 

U-1 

2.83 

11.4 

12.19 

4.16 

35.8 

U-2 

1.63 

10.0 

19.11 

3.25 

53.5 

U-2 

3.0 

162.8 

19.11 

3.25 

7.88 

V-1 

1 .63 

24.6 

13.46 

4.56 

40.1 

V-4 

1.63 

72.3 

14.84 

3.80 

38.8 

H-1 

3.0 

138.8 

29.14 

2.03 

8.43 

S-1 

1.57 

18.4 

58.52 

4.65 

8.97 

S-1 

1 .85 

25.1 

58.52 

4. .65 

8.43 

S-1 

2.1 

32.4 

58.52 

4.65 

7.73 

S-2 

1 .57 

21.3 

89.17 

3.38 

8.04 

S-2 

1.85 

31 .4 

89.17 

3.38 

6.98 

S-3 

1.85 

53.3 

87.99 

3.25 

4.31 

S-3 

2.1 

69.4 

87.99 

3.25 

3.93 

S-4 

1 .57 

49.8 

86.24 

3.13 

3.85 

S-4 

1 .85 

41.1 

86.24 

3.13 

5.95 

S-4 

2.1 

51.4 

86.24 

3.13 

5.64 

S-5 

2.1 

43.3 

60.08 

5.12 

5.14 

S-6 

2.1 

15.6 

59.98 

5.12 

4.61 

S-7 

1.57 

16.8 

58.69 

4.64 

3.16 

S-8 

2.1 

45.9 

93.24 

4.10 

4.46 

S-9 

2.1 

6.4 

92.18 

4.11 

10.47 

S-10 

1 .57 

32.0 

117.65 

4.12 

2.11 


Notes : 1 . 

GP not corrected for C 

^/ 1<2 term. 



2. 

GP corrected 

for C^/Kg 

term. 




TABLE D-2 


PANEL FLUTTER DATA FROM REFERENCE 18 
Isotropic panel supported at all edges 


Data shown in figure 13 

E = 69 GN/m^ 
h =0.81 mm 
a = 0.762 m 
b = 0.17 m 
GP = 4.48 


Panel 

No. 

Log 

No. 

M 

q , 

(kPa) 

FP 

(xlO-4) 

10 

17 

1 .4 

30.6 

2.9 

10 

19 

1 .3 

32.3 

2.8 

10 

20 

1.2 

33.5 

3.1 

10 

21 

1 .2 

38.3 

2.7 

10 

23 

1.2 

38.3 

2.7 

10 

24 

1.1 

55.1 

2.7 

4 

19 

1 .4 

30.4 

2.9 

4 

21 

1 .3 

31 .1 

2.9 

4 

26 

1 .2 

40.6 

2.6 

4 

57 

1 .4 

35.9 

2.5 

4 

70 

1.4 

32.3 

2.8 

5 

6 

1 .4 

33.0 

2.7 

5 

7 

1.33 

32.0 

2.8 

5 

8 

1 .2 

40.7 

2.5 

5 

9 

1.35 

32.3 

2.8 

5 

16 

1 .3 

30.9 

2.9 

5 

18 

1 .4 

32.3 

2.8 

5 

19 

1 .3 

29.9 

3.0 

5 

20 

1 .3 

34.5 

2.6 

5 

21 

1 .2 

35.9 

2.9 

5 

25 

1 .3 

29.4 

3.0 

5 

26 

1 .3 

32.3 

2.8 

5 

27 

1 .4 

31 .1 

2.9 

5 

28 

1 .4 

34.5 

2.7 

5 

36 

1 .3 

32.1 

2.8 

5 

37 

1 .2 

38.1 

2.7 

5 

38 

1 .3 

32.3 

2.8 

5 

54 

1 .3 

35.4 

2.6 

5 

82 

1 .3 

37.8 

2.4 

5 

83 

1 .2 

43.8 

2.4 

5 

84 

1 .3 

39.7 

2.3 
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TABLE D-2 (Concluded) 

PANEL FLUTTER DATA FROM REFERENCE 18 


Panel 

No. 


Log 

No. 

M 

q 

(kPa) 

(xlO-4) 

6 

1.4 

35.2 

2.6 

7 

1.3 

34.2 

2.6 

8 

1.2 

40.0 

2.7 

14 

1.3 

39.5 

2.3 

35 

1 .3 

41.2 

2.2 

36 

1.3 

41 .4 

2.2 
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TABLE D-3 

PANEL FLUTTER DATA FROM REFERENCE 21 


Note : 


Isotropic panel supported at all edges 
Data shown in figure 13 


a = 0.216 m 
D = 0.75 N-m 


Panel 10-20-20 


b 

GP 


0.470 m 
0.46 


M 

2.0 

3.0 


q FP 

(kPa) (x 10"3) 

29.44-36.82 4.40-3.52 

59.29-61.71 3.57-3.43 


Panel 10-20-16 


b 

GP 


a = 0.216 m 
D = 0.39 N-m 

M 


2.0 

3.0 


q 

(kPa) 

17.17-19.65 

23.72-29.65 


= 0.470 m 
= 0.46 

FP 

(x 10-3) 

3.87-3.37 

4.57-3.66 


a = 0.216 m 
D = 0.16 N-m 


Panel 10-20-12 


b 

GP 


= 0.470 m 
= 0.46 


M 

q 

(kPa) 

2.0 

7.38- 8 

3.0 

11.86-14 

4.0 

20.14-22 

5.0 

26.21-27 


FP 

(x 10-3) 

.58 3.80-3.26 

.83 3.79-3.03 

.62 3.05-2.72 

.31 2.97-2.85 


The q data presented in this table were taken from figure 16 
of reference 21 . 
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TABLE D-3 (Continued) 

PANEL FLUTTER DATA FROM REFERENCE 21 


Panel 15-20-20 
a = 0.343 in b 

D = 0.75 N-m GP 


= 0.470 m 
= 0.73 


M 


2.0 

2.5 

3.0 


q 

(kPa) 


(x 10'^) 


11.03-12.27 2.93-2.64 

14.20- 17.79 3.01-2.41 

14.20- 16.62 3.72-3.18 


Panel 20-15-20 
a = 0.470 m b 
D = 0.75 N-m GP 


= 0.343 m 
= 1.37 


M 


2.0 

3.0 

4.0 
4.5 


q FP 

(kPa) (x 10-3) 


8.58- 9.79 
10.68-11 .86 
17.58-20.13 
15.51-17.24 


1 .47-1 .28 
1 .92-1 .73 
1 .60-1 .40 
2.05-1 .85 


Panel 20-10-25 
a = 0.470 m b 
D = 1 .47 N-m GP 


= 0.216 m 
= 2.18 


M 


2.0 

2.5 

3.0 


q 

(kPa) 

9.79-11 .03 
17.79-22.20 
23.72-26.61 


FP 

(x 10-3) 

2.51-2.23 
1.83-1.46 
1 .69-1 .51 


Panel 20-10-20 
a = 0.470 m b 
D = 0.75 N-m GP 


= 0.216 m 
= 2.18 


M 


2.0 

3.0 


q 

(kPa) 

14.75-19.65 

23.72-26.68 


(x 10"4) 


8.51-6.40 

8.66-7.70 
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TABLE D-3 (Concluded) 

PANEL FLUTTER DATA FROM REFERENCE 21 


a = 0.470 m 
D = 0.39 N-m 


Panel 20-10-16 


b = 0.216 m 
GP = 2.18 


M 

q 



(kPa) 

(x 10-4) 

2.0 

7.38- 9.79 

8.73- 6.58 

2.5 

8.89-10.69 

9.58- 7.97 

3.0 

9.48-11.86 

11.09- 8.87 

3.5 

11.31-13.17 

11.03- 9.47 

4.0 

13.86-15.10 

10.39- 9.53 

4.5 

11.24-12.06 

14.51-13.52 

5.0 

22.75-25.02 

8.00- 7.28 


Panel 20-06-12 


a = 0.470 m 


b = 0.114 m 

D = 0.16 N-m 


GP = 4.11 

M 

q 

FP 


(kPa) 

(x 10-4) 

2.0 

12.27-14.75 

2.21-1 .84 

3.0 

20.75-23.72 

2.14-1.87 

4.0 

25.17-30.20 

2.41-2.01 
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TABLE D-4 




Note: 


PANEL FLUTTER DATA FROM REFERENCE 8 
Orthotropic panel with flow angularity 
Data shown in figure 15 



a = 0. 

48 m 



b = 0. 

48 m 



Klt = 1 

.94 



D12/D1 

= 0.17 



Dg/Di 

= 0.00026 

A 

X 

GP 

FP 

(deg) 



(x 10-4) 

30 

11 .50 

3.26 

1 .42 

45 

8.90 

3.32 

1 .12 

60 

7.05 

3.34 

1 .02 

75 

7.05 

3.34 

.86 

90 

8.16 

3.35 

.64 

The above 

data were 

taken 

from figure 6 

reference 

8 wi th a 

normal 

ization factor 

on X of 371 where 




_ 2qa^ 

^ ■ D^?W 
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Panel 


Note : 


TABLE D-5 

PANEL FLUHER DATA FROM REFERENCE 11 
Orthotropic panel with flow angularity 

Data shown in figures 15 and 16 


a 

(m) 

a/b 

‘^2 

(N-m) 

D 1 /D 2 

D 12 /D, 

2 \t 

0.572 

5. 

0.0982 

755 

1.5 

Clamped 

.356 

.933 

.405 

420 

137. 

1.13 

.356 

.933 

.405 

420 

137. 

1.13 

.356 

.933 

.405 

420 

137. 

.252 

.356 

.933 

.405 

270 

137. 

.698 




Panel 1 



A 


q 

GP 


FP 

(deg) 


(kPa) 



(xlO'^) 

31.6 


77.0 

0.23 


3.8 

32.5 


72.7 

.23 


3.9 

33.5 


68.4 

.23 


4.0 

36.2 


64.2 

.23 


3.9 

46. 


56.4 

.23 


3.4 

46. 


60.0 

.24 


3.2 

50. 


58.1 

.24 


3.0 

57. 


53.0 

.24 


2.9 


The q's presented in this table were calculated from q/f(M) data 
provided by Peter Shyprykevich, Grumman Aerospace Corporation. 
This information is summarized in figures 11 through 15, 
reference 11 . 
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TABLE D-5 (Continued) 

PANEL FLUTTER DATA FROM REFERENCE 11 


Panel 2 


A 

(deg) 

(kPa) 

GP 

FP 

(XIO"^) 

19. 

61.5 

1.81 

1.13 

19.8 

68.3 

1.82 

.97 

19.8 

56,4 

1.82 

1.17 

21.7 

51 .2 

1.84 

1.16 

24.2 

42.6 

1,87 

1.23 

28.8 

34.2 

1,90 

1.24 

30.2 

25.6 

1,91 

1 .57 

32. 

25.6 

1.92 

1 .47 

33.1 

34.2 

1.92 

1 .05 

33.2 

22.2 

1.92 

1.62 

37.2 

17.0 

1.93 

1 .84 

59.5 

10.3 

1.96 

1 .71 


Panel 3 


A 

(deg) 

q 

(kPa) 

GP 

FP 

(xlO""^) 

35.6 

85.5 

1.93 

3.87 

36.1 

76.7 

1.93 

4.23 

36.5 

85.5 

1.93 

3.75 

39. 

68.2 

1.94 

4.33 

41.8 

59.7 

1.94 

4.53 

46.8 

59.7 

1.95 

3.94 

49.1 

51.3 

1.95 

4.33 

54. 

42.8 

1.96 

4.61 

64.8 

34.3 

1.97 

4.66 

69.4 

30.7 

1*97 

4.88 

90.5 

27.4 

1,97 

4.74 


Panel 4 


A 

(deg) 

q 

(kPa) 

GP 

FP 

(xlO"^) 

16.3 

25.7 

0.92 

4.27 

17.1 

34.3 

.92 

3.02 

17.9 

17.1 

.92 

5.73 

18.5 

42.8 

.92 

2.20 

19.3 

51 .4 

.93 

1 .74 

21.2 

17.1 

.93 

4.64 

25.4 

8.6 

.93 

7.47 
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TABLE D-5 (Concluded) 

PANEL FLUTTER DATA FROM REFERENCE 11 


Panel 5 


A 

q 

GP 

FP 

(deg) 

(kPa) 


(xio"^) 

15.5 

85.6 

1.17 

1.11 

17. 

77.0 

1.18 

1.11 

20.5 

68.4 

1.19 

1.01 

23. 

60.0 

1.19 

1.01 

26. 

51.4 

1.20 

1.02 

30. 

42.8 

1.20 

1.03 

40.5 

34.3 

1.21 

.83 

58.3 

25.7 

1.21 

.68 



I 


TABLE D-6 

PANEL FLUTTER DATA FROM REFERENCE 12 



Isotropic panel 

supported at 

all 

edges 


with 

inplane loads 



Data shown in f' 

igures 

17 and 

21 




Panel 

1 




a = 0.343 

m 

b = 0. 

343 

m 


D = 7.085 

N-m 

M = 1. 

96 




Kj = 968 



^CR 

, ^ V 

(kPa) 


GP 


FP 

(xlO“^) 

0.13 

66.6 


0.932 


2.66 

.17 

61.1 


.911 


2.26 

.18 

61.7 


.905 


2.11 

.24 

51.4 


.871 


1 .76 

.27 

51 .4 


.854 


1.49 

.48 

33.7 


.721 


.87 

.50 

33.7 


.707 


.81 

.57 

20.0 


.655 


1 .07 

.58 

20.6 


.648 


1 .01 

.70 

14.6 


.547 


1 .00 

.71 

14.6 


.538 


.97 

.92 

8.6 


.283 


1 .00 


Panel 2 
a = 0.622 m 
D = 13.839 N-tn 


*^08 

q 

(kPa) 

0.03 

42.5 

.29 

22.6 

.56 

15.6 

.63 

8.8 

.71 

6.2 

.95 

6.2 


b = 0.318 m 


M = 1.96 


43.4 


GP 

FP 


(xlO"^) 

1 .90 

1.77 

1 .63 

1.10 

1.28 

.93 

1 .18 

1 .48 

1.04 

1 .89 

.43 

1 .45 


Note: For all panels the Pj,|^ and q data were obtained from information shown 

in figure 9, reference 12. 
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TABLE D-6 (Continued) 

PANEL FLUTTER DATA FROM REFERENCE 12 





Panel 

3 




a = 0. 

,622 

m 

b = 0, 

.257 m 



D = 13.839 N-m 

M = 1, 

.96 






22.8 



p 

CR 

q 

(kPa) 



GP 


FP 

(xlO"^) 

0.21 

65.7 



2.11 


0.52 

.46 

47.4 



1.74 


.45 

.53 

30.4 



1.62 


.64 

JZ 

18.6 



1.25 


.85 

.90 

11.9 



.75 


1 .14 

.97 

7.6 



.41 


1.68 




Panel 

4 




a = 0. 

622 

m 

b = 0. 

.216 m 



0 = 7. 

085 

N-m 

M = 1, 

.96 






25.7 



Ppn 

q 



6P 


FP 

CK 

(kPa) 





(xlO"^ 

0.13 

51.6 



2.64 


0.45 

.15 

51.6 



2.60 


.42 

.27 

33.6 



2.41 


.51 

.29 

33.6 



2.38 


.49 

.47 

22.3 



2.06 


.58 

.56 

17.8 



1.87 


.67 

.88 

9.2 



.98 


1.01 

.99 

9.1 



.28 


.95 
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TABLE D-6 (Continued) 

PANEL FLUTTER DATA FROM REFERENCE 12 


Panel 5 



a 

= 0.622 in 

b = 0.165 m 



D 

= 3.628 N-m 

M = 1.96 




S = 

21 .6 


p 

* no 


q 

GP 

FP 

CK 


(kPa) 


(x10“^) 

0.25 


68.6 

3.19 

1.59 

.38 


60.4 

2.90 

1.58 

.38 


51 .2 

2.90 

1.86 

.39 


60.4 

2.88 

1.57 

.41 


34.1 

2.83 

2.73 

.43 


50.4 

2.78 

1 .82 

.50 


23.8 

2.61 

3.65 

.52 


23.9 

2.55 

3.59 

.68 


17.0 

2.08 

4.58 

.70 


17.1 

2.02 

4.51 

.88 


8.8 

1.28 

8.03 



Panel 6 



a 

= 0.622 m 

b = 0.148 m 



D 

= 3.628 N-m 

M = 1.96 





14.3 


P pp 


q 

GP 

FP 

LK 


(kPa) 


(xlO“^) 

0.10 


60.3 

3.85 

2.41 

.38 


43.1 

3.20 

2.44 

.41 


51 .7 

3.12 

1 .99 

.49 


34.6 

2.90 

2.84 

.64 


26.3 

2.44 

3.44 

.66 


26.3 

2.37 

3.41 

.68 


17.7 

2.30 

5.02 

.69 


18.0 

2.26 

4.91 

.78 


18.0 

1 .91 

4.73 

.94 


14.3 

1 .00 

5.61 

.95 


14.3 

.91 

5.59 
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TABLE D-6 (Continued) 
PANEL FLUTTER DATA FROM REFERENCE 12 

Panel 7 



a = 0.343 

m 

b = 0.343 m 



D = 13.839 

N-m 

M = 1.96 





545 


Pro 

q 


GP 

FP 

CR 

(kPa) 



X 

o 

I 

CO 

0.16 

57.9 


0.916 

4.97 

.39 

63.8 


.780 

1.29 

.39 

57.9 


.780 

1.43 

.42 

52.0 


.761 

1.40 

.51 

34.1 


.699 

1.50 

.52 

17.5 


.692 

2.83 

.54 

42.6 


.678 

1.08 

.58 

26.0 


.647 

1.56 

.66 

17.5 


.583 

1 .82 

.69 

12.6 


.556 

2.31 

.77 

11.7 


.479 

2.02 

.90 

13.8 


.316 

1.27 



Panel 

8 



a = 0.622 m 

b = 0.318 m 



D = 29.007 

N-m 

M = 1.96 





21.6 


Ppn 

q 


GP 

FP 

CK 

(kPa) 



(xlO"^) 

0.23 

75.8 


1 .68 

0.81 

.24 

42.8 


1.67 

1.40 

.26 

42.8 


1.64 

1.32 

.27 

58.0 


1.63 

.94 

.29 

57.7 


1.61 

.90 

.63 

22.6 


1.16 

1.21 

.65 

22.6 


1.13 

1 .17 

.70 

17.5 


1.05 

1.42 

.72 

22.0 


1.01 

1.10 

.90 

8.6 


.60 

2.30 

.92 

8.6 


.54 

2.26 
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TABLE D-6 (Continued) 



PANEL FLUTTER DATA 

FROM REFERENCE 

12 


Panel 

9 



a = 0.622 m 

b = 0.216 m 



D = 13.839 N-m 

M = 1.96 



= 10.2 


p 

q 

GP 

FP 

CR 

(kPa) 


(xlO“^) 

0.09 

76.5 

2.62 

0.67 

.12 

76.2 

2.58 

.61 

.15 

76.2 

2.53 

.56 

.39 

60.7 

2.15 

.46 

.40 

60.1 

2.13 

.46 

.44 

40.4 

2.06 

.65 

.46 

40.4 

2.02 

.64 

.53 

40.7 

1.88 

.59 

.58 

20.2 

1.78 

1 .13 

.60 

29.9 

1.74 

.75 

.63 

19.9 

1 .67 

1 .09 

.67 

19.9 

1.58 

1 .06 

.69 

29.9 

1 .53 

.69 

.71 

8.5 

1.48 

2.40 

.73 

20.2 

1 .43 

.99 

.74 

19.9 

1 .40 

1.00 

.76 

12.6 

1.35 

T.56 

.78 

12.6 

1 .29 

1 .54 

.81 

7.9 

1 .20 

2.40 


Panel 10 

a = 0,622 tn b = 0.165 m 

D = 7.085 N-m M = 1.96 

K; - 8.45 


p 

^^CR 

q 

(kPa) 

GP 

FP 

(xlO'2) 

0.05 

51.2 

3.47 

0.62 

.14 

43.1 

3.31 

.58 

.37 

35.0 

2.83 

.54 

.53 

25.5 

2.44 

.65 

.55 

25.2 

2.39 

.65 

.58 

17.7 

2.31 

.91 

.61 

25.2 

2.23 

.63 

.61 

17.7 

2.23 

.89 

.67 

9.0 

2.05 

1 .70 

.67 

12.2 

2.05 

1 .25 

.68 

25.4 

2.02 

.60 

.78 

11.6 

1 .67 

1 .25 
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TABLE D -6 (Continued) 

PANEL FLUTTER DATA FROM REFERENCE 12 


Panel 11 

a = 0.662 m b = 0.216 m 

D = 7.085 N^m M = 1.96 

K 3 = 22.3 


p 

q 

GP 

FP 

CR 

(kPa) 


(xlO’^) 

0.01 

35.0 

2.98 

0.97 

.02 

30.1 

2.97 

1.03 

.05 

30.1 

2.92 

.85 

.13 

20.6 

2.80 

.95 

.27 

13.6 

2.56 

1.09 

.42 

7.9 

2.28 

1.54 

.43 

9.6 

2.26 

1.25 

.43 

7.7 

2.26 

1.56 

.49 

7.7 

2.14 

1.47 

.62 

7.0 

1.85 

1.44 

.74 

8.2 

1.53 

1.13 

.86 

8.5 

1.12 

1.00 

.86 

9.5 

1.12 

.90 



Panel 12 




a = 0.622 m b = 

0.165 m 



D = 7.085 N-m M = 

1.96 



Kg = 9.96 



P 

'^CR 

q 

GP 

FP 

(kPa) 


(xlO'^) 

0.06 

82.3 

3.48 

3.71 

.08 

82.3 

3.45 

3.49 

.13 

65.3 

3.35 

3.92 

.15 

69.3 

3.31 

3.55 

.17 

68.7 

3.27 

3.47 

.49 

52.6 

2.57 

3.25 

.49 

65.1 

2.57 

2.63 

.51 

54.2 

2.52 

3.11 

.53 

54.2 

2.46 

3.07 

.53 

61.9 

2.46 

2.69 

.56 

52.7 

2.38 

3.10 

.57 

53.8 

2.36 

3.01 

.58 

34.4 

2.33 

4.68 
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TABLE D-6 (Concluded) 

PANEL FLUTTER DATA FROM REFERENCE 12 


Panel 12 (Concluded) 


p 

*^CR 

(k?a) 

GP 

FP 

(xlO""^) 

.58 

53.8 

2.33 

2.99 

.60 

34.3 

2.27 

4.64 

.72 

26.5 

1.90 

5.62 

.76 

26.2 

1.76 

5.57 

.78 

16.5 

1 .68 

8.76 

.80 

17.0 

1.61 

8.42 

.84 

18.4 

1 .44 

7.64 

.85 

18.7 

1.39 

7.48 




Panel 13 



a = 0.622 

m b = 0.216 m 



D = 7.085 

N-m M = 1.96 




Ks = 14.1 


p 

•^CR 

q 

(kPa) 

GP 

FP 

(xlO"^) 

0.11 

41.6 

2.63 

0.59 

.17 

34.4 

2.54 

.60 

.21 

38.4 

2.47 

.50 

.30 

26.3 

2.33 

.62 

.37 

17.8 

2.21 

.82 

.54 

8.8 

1 .89 

1.37 
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Note: 


Note: 


TABLE D-7 

PANEL FLUTTER DATA FROM REFERENCE 7 

Orthotropic panels with spring supported stream 
edges and inplane loads 

Data shown in figure 18 

Panel I 


a = 0.605 m b = 0.605 n 

D, = 1 N-m D, = 1640 N-m 

D ' = 376 N-m Kc = 63 

It *1 ^ 



12 

M = 3 



Test 

Pro 

q 

GP 

FP 


CR 

(kPa) 


(xlO-^^) 

1 

0.74 

no 

6.42 

1.15 

1 

.92 

101 

3.56 

1.24 

2 

.40 

134 

9.76 

.95 

2 

.71 

115 

6.79 

1.10 

4 

.72 

118 

6.67 

1.07 

5 

.43 

137 

9.51 

.93 

6 

.61 

106 

7.87 

1.19 

8 

.80 

107 

5.63 

1.18 

determi ned 

from 

information 

shown in 

figure 8 , reference 7. 


Panel II 

a = 0.605 m b = 0.605 m 

D, = 0.94 N-m D. = 4150 

= 798 N-m K 3 = 6.9 


Test 

^CR 

q 

(kPa) 

GP 

FP 

(xlO’"^) 

15 

0.38 

98 

4,62 

1.22 

16 

.09 

131 

5.60 

.92 

17 

.45 

117 

4.35 

1.02 

17 

.69 

102 

3.27 

1.17 

18 

.86 

83 

2.20 

1.44 

determined from 

information 

shown in 

figure 9, reference 7 
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TABLE D-8 

PANEL FLUTTER DATA FROM REFERENCE 19 


Isotropic panel with inplane loads 

Data shown in figures 18 and 20 

E = 72.4 GN/m^ 
a = 0,66 m 
M = 3 


Partially Clamped Panel 
a/b = 3.3 


h 

p 

CR 

q 

GP 

FP 

(mm) 

(kPa) 


(xlO'^) 

1.4 

0.44 

117 

2.47 

0.43 

1.4 

.62 

69 

2.03 

.63 

2.03 

.63 

158 

2.01 

.84 

2.03 

.78 

138 

1.55 

.88 

2.03 

.77 

117 

1.58 

1.04 

2.03 

.85 

83 

1.28 

1.40 

2.03 

.94 

69 

.81 

1.61 



Partially Clamped 

Panel 




a/b = 3.3 



h 

p 

“^CR 

q 

GP 

FP 

(mm) 

(kPa) 


(xlO"^) 

1.32 

0.39 

179 

2.58 

0.25 

1.37 

.56 

117 

2.19 

.37 

1.37 

.29 

234 

2.78 

.24 

1.37 

.39 

165 

2.58 

.30 

1.40 

.29 

234 

2.78 

.26 

1.60 

.59 

186 

2.11 

.36 

2.59 

.85 

220 

1.28 

1.09 

2.59 

.88 

165 

1.14 

1.44 

2.59 

.93 

117 

.87 

1.99 

2.59 

.95 

96 

.74 

2.38 

2.59 

.96 

69 

.66 

3.30 


Note: based on information shown in figures 10 and 11, reference 19. 
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TABLE D-8 (Concluded) 



PANEL 1 

FLUTTER DATA 

FROM REFERENCE 19 




Fully Clamped Panel 




a/b = 

= 3.3 


h 

Ppn 

q 

GP 

FP , 

(mm) 

CR 

(kPa) 


(xlO’^) 

1.35 

0.26 

186 

2.84 

0.30 

1.35 

.46 

no 

2.42 

.40 

1.35 

.56 

69 

2.19 

.59 

1.37 

.61 

83 

2.06 

.50 

1.60 

.76 

69 

1.62 

.87 

1.60 

.69 

96 

1.84 

.65 

1.65 

.70 

no 

1.81 

.62 

1.65 

.47 

165 

2.40 

.49 

1.93 

.74 

no 

1.68 

.96 

1.93 

.81 

96 

1.44 

1.06 

2.59 

.79 

172 

1.51 

1.45 



Fully Clamped Panel 




a/b = 

3.7 


h 

p 

q 

GP 

FP . 

(mm) 

CR 

(kPa) 


(xlO'^) 

1.37 

0.36 

269 

2.96 

0.21 

1.37 

.55 

200 

2.48 

.25 

1.37 

.65 

138 

2.19 

.34 

1.37 

.71 

117 

1.99 

.39 

1.62 

.57 

269 

2.43 

.30 

1.62 

.78 

96 

1.74 

.75 

1.62 

.74 

83 

1.89 

.89 

1.62 

.87 

69 

1.34 

1.01 

1.62 

.88 

76 

1.28 

.91 
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TABLE D-9 

PANEL FLUTTER DATA FROM REFERENCE 23 
Supported isotropic panel with inplane loads 
Data shown in figure 18 


M = 3 

E = 72.4 GN/m^ 


Panel 

Test 

h 

a 

b 

P 

q 

GP 




(mm) 

(m) 

(m) 

CK 

(kPa) 


(XIO"^) 

2-A 

1 

0.64 

0.686 

0.0686 

0.58 

177 

6.48 

5.78 

2-A 

5 

.81 

.686 

.0686 

.86 

212 

3.74 

9.59 

2-A 

6 

.81 

.686 

.0686 

.86 

216 

3.74 

9.41 

4-A 

8 

.64 

.66 

.066 

.46 

240 

7.35 

4.84 

4-A 

9 

.64 

.66 

.066 

.71 

144 

5.39 

7.89 





M = 3 








E 

= 200 GN/rr 

,2 




Panel 

Test 

h 

a 

b 

PpD 

q 

GP 

PP t; 



(mm) 

(m) 

(m) 

CR 

(kPa) 


(xlO"^) 

4-S 

3 

0.51 

0.711 

0.0711 

0.68 

237 

5.66 

5.37 

4-S^ 

4 

.51 

.711 

.0711 

.69 

237 

5.57 

5.37 

4-Sp 

5 

.51 

.711 

.0711 

.66 

199 

5.83 

6.41 

4 -S 2 

6 

.51 

.711 

.0711 

.15 

231 

9.23 

5.93 


Note: Pqi^'s based on information presented in table I, reference 23. 
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TABLE D-10 

PANEL FLUTTER DATA FROM REFERENCE 24 
Clamped isotropic panel with inplane loads 
Data shown in figure 18 


a = 0.635 m 
b = 0.159 m , 
E = 72.4 GN/m'^ 
M = 3.0 


Test 

h 

(mm) 

p 

*^CR 

q 

(kPa) 

GP 

3 

(XIO 

2 

1.6 

0.20 

240.8 

3.58 

0.54 

3 

1.6 

.95 

166.1 

.89 

.48 

4 

1.6 

.95 

119.2 

.89 

.67 

5 

2.06 

.51 

214.0 

2.80 

.99 

6 

2.06 

.88 

78.5 

1.39 

2.24 

7 

3.18 

1.00 

236.5 

0.0 

2.61 

Determination 
reference 24. 


based on 

information given 

in figure 

14, 
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TABLE D-n 

PANEL FLUTTER DATA FROM REFERENCE 25 


Note: 


Note: 


Supported isotropic panel with inplane loads 
Data shown in figure 18 




a = 0.661 m . 






E = 72.4 GN/m' 
M = 3 






a/b = 2.5 




Test 

h 

^CR 

q 

GP 

FP o 


(mm) 

(kPa) 


(xlO"'^ 

1 

1.6 

0.96 

129.5 

0.50 

0.26 

2 

1.6 

.56 

187.2 

1.66 

.26 

3 

1.8 

.29 

141.0 

2.11 

.73 

4 

1.8 

.39 

93.8 

1.95 

.92 

5 

1.8 

.52 

70.9 

1 .73 

1.03 

6 

2.62 

.56 

151.3 

1.66 

1.41 

8 

2.62 

.65 

140.3 

1.48 

1.38 

10 

2.62 

.77 

232.0 

1.20 

.75 

11 

2.62 

.65 

174.3 

1.48 

1.11 

CR ^ 

based on information 

presented on 

figure 7, 

reference 

25. 



a/b = 2.9 




Test 

h 

P 

^CR 

q 

GP 



(mm) 

(kPa) 


(xlO"^ 

1 

2.06 

0.52 

235.6 

2.01 

5.49 

2 

2.06 

.77 

187. 

1.39 

5.61 

3 

2.06 

.88 

128.1 

1.00 

7.61 

6 

2.06 

.98 

111.8 

.41 

8.20 

7 

1.35 

.50 

117.5 

2.05 

3.15 

8 

1 .35 

.77 

80.4 

1.39 

3.66 

9 

1.35 

.09 

152.7 

2.77 

5.02 

10 

1.6 

.62 

146.3 

1.79 

3.78 

11 

1.6 

.31 

146. 

2.41 

5.36 

12 

1.6 

.69 

81.9 

1.61 

6.37 

13 

1.6 

.34 

176. 

2.36 

4.25 


^Cr'^ based on information presented on figure 8, reference 25, 
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TABLE D-12 

PANEL FLUTTER DATA FROM REFERENCE 26 


Orthotropic panels with spring supported stream 
edges and inplane loads 

Data shown in figure 18 

Panel III 

a = 0.605 m b = 0.605 

D, = 2.4 N-m Do = 2130 N-m 

D,i = 500 N-m IC = 89 

M = 3 ^ 


Test 

p 

•^CR 

q 

(kPa) 

6P 

FP 4 
(xlO-^) 

2 

0.54 

204 

7.86 

1.48 

3 

.42 

221 

8.83 

1.38 

4 

.44 

225 

8.68 

1.35 

5 

.36 

220 

9.28 

1.39 

6 

.42 

213 

8.83 

1.43 

7 

.30 

232 

9.70 

1.32 

7 

.78 

179 

5.44 

1.67 

7 

.78 

201 

5.44 

1.49 

7 

.75 

212 

5.79 

1.41 


Note: Pqp's 


determined from information 


shown in figure 8(a), reference 26 


Panel 

a = 0.483 m 
D, = 0.87 N-m 
D,p = 556 N-m 

M = 


b = 0.483 
Do = 3380 
1C = 2.1 


Test 

^CR 

q 

(kPa) 

GP 

FP 

(xlO 

1 

0.77 

154 

1.70 

1.41 

1 

.85 

141 

1.37 

1.54 

2 

.49 

190 

2.53 

1.15 

3 

.68 

153 

2.00 

1.42 

4 

.50 

213 

2.50 

1.02 

5 

.39 

191 

2.76 

1.14 

5 

.96 

95 

.71 

2.28 


Note: P(^(^'s determined from information shown in figure 8(b), reference 26 
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TABLE D-13 


PANEL FLUTTER DATA FROM REFERENCE 10 


Orthotropic panel with stream edge support flexibility 
and pressure differential effects 


Data shown in figure ig 


Panel V-5 

a = 0.61 m b = 0.61 m 

D, = 2.4 N-m Dp = I960 N-m 

D ^2 = 433 N-m 1^ = 4.3 

E = 73 GN/m^ 


M 

Ap 

(kPa) 

(kPa) 

GP 

FP^^4 

(xlO-^) 

(xlO"^) 

1.63 

0.24 

14.2 

4.19 

8.81 

1.26 

1.84 

-.24 

22.1 

4.19 

7.27 

1.04 

1.96 

-.29 

24.5 

4.19 

7.25 

1.16 




Panel V-2 




a 

= 0.61 m b = 0.61 

m 



°i 

^12 

=1.46 N-m Dp = 3560 

= 464 N-m = 1.2 

N-m 




E = 214 GN/m"^ 



M 

Ap 

q GP 

FP^^J 

Fp(2) 


(kPa) 

(kPa) 

(xlO'4) 

(xlO-3) 

1 .63 

-0.24 

27.80 2.99 

2.73 

0.697 

1.84 

.24 

28.01 2.99 

3.49 

.890 

1.96 

- .24 

24.08 2.99 

4.49 

1.144 


Notes: 1. FP not corrected for Ap effects. 
2. FP corrected for Ap effects. 
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TABLE D-14 


PANEL FLUTTER DATA FROM REFERENCE 27 


Isotropic panel at low supersonic flow 
Data shown in figure 22 


M 

a = 0.229 m 
b = 0.457 m 
h = 1 .019 mm p 
E = 38.6 GN/m^ 

6 

q 


(mm) 

(kPa) 

1.05 

24.13 

50.3 

1.1 

26.16 

49.6 

1.1 

7.62 

29.2 

1.2 

25.9V 

48.7 

1.2 

19.56 

42.3 

1.2 

13.21 

36.7 

1.2 

11.68 

34.2 

1.2 

8.64 

32.9 

1.2 

7.11 

29.0 

1.3 

24.38 

52.6 

1.3 

7.11 

42.2 

1.4 

22.61 

61.8 

1.4 

6.86 

56.9 


Note: The quantity 6 is the boundary layer thickness at the test conditions 



I 


TABLE D-15 

PANEL FLUTTER DATA FROM REFERENCE 28 


Isotropic panel at low supersonic flow 
Data shown in figure 22 



a = 0.229 m 



b = 0.457 m 



h = 0.488 mm o 



E = 144.8 GN/m^ 


M 

6 

q 


(mm) 

(kPa) 

1.2 

28.96 

19.3 

1.2 

20.83 

15.6 

1.2 

16.76 

12.3 

1.2 

11.68 

10.6 

1.3 

29.21 

22.2 

1.3 

9.65 

13.4 

1.3 

8.13 

12.7 

1.4 

11.94 

18.7 

1.4 

8.64 

17.7 


Note: The quantity 6 is the boundary layer thickness at the test conditions. 
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